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By W. G. Lidman and A. R. Bobrovsky 


SUMMARY 

An analysis is made to determine which properties of materials 
affect their resistance to fracture by thermal stresses. From this 
analysis, a parameter is evaluated that is correlated with the 
resistance of ceramic materials to fracture by thermal shock as 
experimentally determined. This parameter may be used to predict 
qualitatively the resistance of a material to fracture by thermal 
shock. 

Resistance to. fracture by thermal shock is shown to be dependent 
upon the following material properties: thermal conductivity, 

tensile strength, thermal expansion, and ductility modulus. For 
qualitative prediction of resistance of materials to fracture by 
thermal shock, the parameter may be expressed as the product of 
thermal conductivity and tensile strength divided by the product of 
linear coefficient of thermal expansion and ductility modulus of 
the specimen. 


INTRODUCTION 

The operating temperatures of aircraft propulsion systems are 
limited by the melting temperature of the alloys currently avail- 
able. In general, higher operating temperatures result in higher 
cycle efficiencies for gas turbines. One means of obtaining these 
higher efficiencies is to select materials that are more refractory 
than the alloys currently used. Oxides and carbides of metals are 
therefore being considered for use in aircraft propulsion systems 
because of their high melting temperatures. 

In the selection of a ceramic for use as a turbine-blade mate- 
rial, the properties of the material that affect its operating 
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characteristics must "be considered. Because sudden engine starts 
and stops are relatively frequent during service operation, the 
resistance of Blade materials to fracture by thermal shock is of 
considerable importance. The stresses induced in a material by the 
thermal shocks of unsteady engine operation are the result of heat 
flow from the center of the blade to the cooler surrounding atmos- 
phere upon stopping, or from the hot gas to the blade upon starting. 
This time-dependent type of heat flow is commonly termed "unsteady- 
state heat flow," and the stresses induced are different from those 
that would be induced by having steady flow from a constant- 
temperature hot zone to a constant -temperature cold zone (steady- 
state heat flow). 

The methods of evaluating resistance to thermal shock have not 
as yet been standardized. A general method of determining thermal- 
shock resistance of a material consists in subjecting a specimen to 
a number of cycles of alternate heating and cooling until fracture 
occurs. (The variations in the procedure are considerable.) 

Other investigators have presented work showing the dependence 
of material properties on resistance to thermal shock. In refer- 
ence 1, a derivation of the theory is presented for the tendency of 
bricks to spall. It is shown that the spalling tendency is pro- 
portional to the coefficient of thermal expansion divided by the 
product of the maximum tensile strain and the square root of the 
diffusivity. 

According to reference 1, spalling is usually defined as a 
fracture of a refractory brick or block resulting from any of the 
following causes: (1) temperature gradients in the brick, (2) com- 

pression in a structure of refractory bricks sufficient to cause 
shear failures, and (3) variation in coefficient of thermal expan- 
sion between surface layer and body of the brick. The A.S.T.M. 
method of classifying brick consists in expressing resistance to 
spalling as weight loss after a specified spalling test. This 
weight loss is the result of flaking of the surface layer, or 
chipping of the corners of the brick. 

Fracture instead of spalling has been observed in investiga- 
tions dealing with the thin ceramic specimens that are of interest 
as gas-turbine -blade materials. For this reason, resistance of 
thin ceramic specimens to fracture by thermal shock is of greater 
interest than their resistance to spalling. 

An analysis was therefore made at the MCA Lewis laboratory to 
determine which properties of a ceramic material affect the resist- 
ance to fracture by thermal stresses. The stresses caused by heat 
flow are considered in the analysis reported. 
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SYMBOLS 

The following symbols are used in this analysis: 
area, (sq in.) 

specific heat at constant pressure, (Btu/(lb) (°F) ) 
size factor 

ductility modulus of specimen, (lb/sq in.) 
ductility modulus of element A, (lb/sq in.) 
ductility modulus of element B, (lb/sq in.) 
force, (lb) 

force in element A, (lb) 
force in element B, (lb) 
diffueiVity, (sq ft/hr) 

thermal conduct 1 vity , (Btu/ (hr ) ( sq f t ) ( °F/in . ) ) 
length, (in.) 

convection heat -transfer coefficient, (Btu/(hr)(sq ft)(°F)) 

stress in specimen, (lb/sq in.) 

stress in element A, (lb/sq in.) 

stress in element B, (lb/sq in.) 

tensile strength of specimen, (lb/sq in.) 

temperature of specimen, (°F) 

temperature of element A, (°F) 

temperature of element B, (°F) 

temperature of gas surrounding specimen, (°F) 

time, (sec) 
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u velocity, (ft/sec) 

x direction of heat flow 

y, z directions normal to direction of heat flow 
a linear coefficient of thermal expansion, (in./in./°F) 

€ emissivity of body 

H viscosity of gas 

p density, (lb/cu in.) 

a Stefan-Boltzmann constant, 0.173 x 10~ 8 (Btu/(sq ft)(hr)(°R 4 )) 

ANALYSIS 

The thermal stresses Induced during a thermal -shock test are 
the result of temperature gradients from the center of the specimen 
to the surface, and not of temperature variations along the surface. 
An examination of the equations for heat flow for specific boundary 
conditions yields information about the dependence of heat flow on 
some independent measurable material properties . 


Material Properties Affecting Thermal Shock 


For the unsteady- state condition, the equation for heat flow 
in a small element of the material is 


dT h ( d 2 T . S 8 ! 

* Hi?' 


(i) 


where, for the body, 

h = ( 1 | dimensional constant) 

\ pc p/ V ) 


For a thin specimen such as the disk used in this investigation, 
the heat flow is essentially linear and equation (l) becomes 
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For the case where the outer surface is cooled suddenly, reference 2 
shows that the temperature at any point in the hody is a homogeneous 

function of h and t. Moreover, ^ and are homogeneous 

functions of the same variables. This result indicates that the 
diffusivity does not influence the value of the maximum gradient in 
the hody, hut only the time at which it occurs. 


The stresses that are induced in a material hy the temperature 
gradients are considered in determining the strength requirements 
for adequate resistance to fracture hy thermal shock. It is assumed 
that two adjacent elements A and B of a material are subjected 
to unequal temperatures, the two elements are at temperatures T A 

and Tg (higher than T A ), respectively (fig. 1), and all proper- 
ties of A and B are independent of temperature. If elements A 
and B were free to expand, the expansion of B would be greater 

than that of A as a result of the higher temperature of B. If 

these elements are restrained by each other, the expansions of the 
elements A and B must be the same to a first approximation. 

The solution of the problem of determining the thermal stresses in 

two elements of a body, as shown in figure 1, is based on two con- 
ditions: (1) the tensile force in element A equals the compres- 

sive force in element B; and (2) the total elongation of A 
equals that of B. 

In order to satisfy condition (l). 


The thermal expansion of element A is 

aZ(T A - T 0 ) 

The thermal expansion of element B is 

“»<T b - T 0 ) 

The extension of element A caused by 
it is 

V 

A E a 


(3) 


the tensile force exerted on 


where 
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The compression of element B caused by the compressive force 
exerted on it is 

V 

A V 

where 



or 


In order to satisfy condition (2), 

al < T A - T 0> + TI - “ J < T B - T 0> - ^ 

a(T A * T 0> + | + ^ * “<T a - T 0 ) + OAT 


where 


(4) 


Then 






= aAT 


Because elements A and B comprise a free body in static 
equilibrium, 

S A = 


Then 


For a unit distance, 



’ oE at 

>A 2 Ax 
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In the limiting case of infinitesimal distance, 

aE® 

S A = Sb * S ‘ “T* ‘ 5 > 

For a given value of S, the resistance to cracking by thermal 
shock is determined by the tensile strength of the material. In 

dT 

order for S to be low, aE— must be low. 

7 dx 

Fluid Properties Affecting Thermal Shock 

The usual equation for heat conduction at the surface in one 
dimension is 

-k |£ + q(T - Tq) + oe( T 4 - T 0 4 ) = 0 (6) 

where 

q(T - Tq) convection heat loss 

oc (T 4 - Tq 4 ) radiation heat loss 

An exact solution for a problem in heat conduction in the 
unsteady state with accurate radiation boundary conditions has not 
yet been found (reference 3). The analysis presented herein con- 
siders only the shock encountered in cooling. Inasmuch as radi- 
ation is disregarded, the results are semi quanti tat ive only. 

For the one -dimensional case, the convection heat loss q(T - Tq) 
is equal to and the temperature gradient normal to the length 

of the specimen is equal to 

i = | (T - Tq) (7) 

The convection heat-transfer coefficient q is usually 
expressed in terms of the following gas properties: 

q = f (u,d,p,n,k,Cp) 

For a given size and shape of specimen at a given gas velocity 
and specific heat, the convection heat-transfer coefficient may be 
considered independent of temperature (reference 4). 
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From equation (7) for a low temperature gradient ~ at the 

ox 

surface, a low convection heat-transfer coefficient and a high 
thermal conductivity are desirable properties of the materials. 

From this reasoning, for maximum resistance to cracking by 
thermal stresses, the desirable properties of a material are high 
tensile strength, low coefficient of thermal expansion, low modulus 
of elasticity at failure, high thermal conductivity, and low 
convection heat-transfer coefficient. 

The ratio of stress to strain at failure is designated the 
ductility modulus. and can be evaluated when the stress and the 
strain at failure are known, as shown in figure 2. 

From the analysis presented, a material parameter that includes 
the thermal conductivity, tensile strength, thermal expansion, and 
the ductility modulus of the materials can be used to predict 
qualitatively the susceptibility of a material to cracking by thermal 
shock. This parameter can be expressed as — . 

aE 

APPARATUS AND PROCEDURE 

The equipment used to determine the thermal- shock resistance 
of ceramic materials is shown in figure 3. This unit consists of 
a furnace for heating the specimen to the evaluation temperature 
and an air-quenching system for rapid cooling of the specimen. A 
specimen holder made from a high-temperature alloy is used to 
support the specimen and to transport it from the furnace to the 
air-quenching stream and back. The specimens used for evaluation 
in this unit were disks 2 inches in diameter and l/4 inch thick. 

The procedure used in the thermal-shock evaluation (refer- 
ence 5) consisted in subjecting the specimen to a number of cycles ♦ 
of alternate heating and quenching until failure occurred. The 
specimen was introduced into the preheated furnace where it was 
held at the evaluation temperature for 10 minutes, after which it 
was transported to the quenching air stream and held for 5 minutes. 
Preliminary tests indicated that in this procedure the center of a 
specimen was heated from room temperature to 1800° F in approxi- 
mately 6 minutes and was cooled from 1800° to 400° F in 1 — minutes 

2 

and from 400° to 85° F in 2 ^ minutes. Each thermal-shock cycle 

was followed immediately by the next cycle at the same temperature. 

The appearance of a crack constituted failure of the specimen. 

When a specimen had withstood 25 thermal-shock cycles of a given 
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©valuation temperature, the test was continued at an evaluation 
temperature 200° F higher for 25 cycles or until failure, whichever 
occurred first; if no failure occurred, this evaluation procedure 
was continued at the next higher temperature. 

The thermal-shock specimens were inspected visually after each 
quenching. When there was doubt as to the presence of a crack, 
the specimen was examined "by meane of radiography. 


RESULTS AND DISCUSSION 

The order of merit of some ceramic bodies evaluated in thermal 
shock was experimentally established according to the method 
described in reference 5. A correlation of thermal-shock parameter 
with experimental results is presented in table I. 

With the exception of the tensile strengths of the materials 
at 1800° F, which were determined at this laboratory according to 
the methods described in reference 5, the physical properties of 
these materials were obtained from literature. In seme instances, 
where experimental data were unavailable for a particular body, 
approximate values were determined from data that were available on 
similar materials. These cases are noted in table I. For example, 
data on the stabilized zirconium dioxide Zr0 2 containing 6-percent 

calcium oxide CaO were unavailable; data for pure Zr0 2 were con- 
sequently used in evaluating the parameter for this material. 

The evaluation of the parameter — for five materials indi- 

aE 

cates that increase in resistance to fracture by thermal shock 

accompanies increase in In table I, the value for titanium 

carbide TiC is 12,560 and for stabilized zirconium dioxide Zr0 2 

plus calcium oxide CaO, which is much poorer in thermal-shock 
resistance, the value is 562. 

An examination of the type of fracture sustained during the 
thermal-shock evaluations shows that all the materials listed failed 
during the cooling part of the cycle with the exception of the sta- 
bilized zirconium dioxide, which failed in shear during the first 
cycle. The shear failure, as observed on stabilized zirconium 
dioxide, results from compression produced during the heating portion 
of the shock cycle. In general, if the stress for fracture in shear 
is less than one-half that necessary for failure in tension, a shear 
failure occurs before tensile failure.. The tensile strength for 
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stabilized zirconium dioxide at 1800° F is 5400 pounds per square 
inch (table I); this value is used in the determination of 

otE 

Actually, because the stabilized zirconium dioxide disks failed in 
shear, the value of s should be an equivalent tensile strength 
(somewhat less than the value of 5400 lb/sq in. shown in table I) 
equal to twice the shear strength. 

If a shape and size factor were included in the parameter, a 
dimensionless number might be obtained that would permit predicting 
the resistance of a ceramic material to fracture by thermal shock 
in any shape or size. 


CONCLUSION 

The resistance of a ceramic material to fracture by unsteady- 
state heat flow (thermal shock) can be correlated with a material 
parameter that is qualitatively determined from the product of 
thermal conductivity and tensile strength divided by the product of 
linear coefficient of thermal expansion and ductility modulus. 
Additional experimental data are required on convection heat-transfer 
coefficients before a quantitative prediction of the failure of 
ceramic materials by thermal shock can be made. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio. 
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Direction of heat flow 


4 — 

E A ductility modulus of element A, lb/sq in* 

Eg ductility modulus of element B, lb/sq in. 

I length, in. 

S A stress in element A, lb/sq in* 

Sg stress in element B, lb/sq in* 

T a temperature of element A, °P 

Tg temperature of element B, °P 

Tq temperature of gas surrounding specimen, °P 

a linear coefficient of thermal expansion of 
specimen, in./in. /°F 

Figure 1. - Elongation of two elements of body as 

ri T 

result of temperature gradient ~ in body. 
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Ultimate strength, lb/sq in. 
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Figure 2 • - Relation between ductility modulus E 
and stress and strain at failure* 
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Figure 3. - Apparatus used for determinating resistance of materials to fracture by 

thermal shock. 
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